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Abstract 

A principled framework to generalize variational perturbation approximations (VPA's) 
formulated within the ambit of the nonadditive statistics of Tsallis statistics, is in- 
troduced. This is accomplished by operating on the terms constituting the pertur- 
bation expansion of the generalized free energy (GFE) with a variational procedure 
formulated using q- deformed calculus. A candidate q- deformed generalized VPA 
(GVPA) is derived with the aid of the Hellmann-Feynman theorem. The general- 
ized Bogoliubov inequality for the approximate GFE are derived for the case of 
canonical probability densities that maximize the Tsallis entropy. Numerical exam- 
ples demonstrating the application of the q-deformed GVPA are presented. The 
qualitative distinctions between the q-deformed GVPA model vis-d-vis prior GVPA 
models are highlighted. 

Key words: Generalized Tsallis statistics, additive duality, variational 
perturbation approximations, q-deformed calculus, Hellman-Feynman theorem, 
generalized Bogliubov inequality. 

PACS: 05.20.-y; 05.70.-a; 05.90.-hm 04.20.Fy; 04.25.Nx 



1 Introduction 



The generalized (nonadditive) statistics of Tsallis' [1,2] has recently been the 
focus of much attention in statistical physics, and allied disciplines. Nonad- 
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ditive statistic^, which generahzes the extensive Boltzmann-Gibbs-Shannon 
(B-G-S) statistics, has much utihty in a wide spectrum of disciphnes ranging 
from complex systems and condensed matter physics to financial mathematics 
[!]• Recent works have extended the scope of Tsallis statistics, by demonstrat- 
ing its efficacy in lossy data compression in communication theory [3] and 
machine learning [4]. 

Variational perturbation approximations (VPA's) [5] are extensively employed 
in quantum mechanics and statistical physics [6]. The first attempt to general- 
ize VPA's to the case of Tsallis statistics was performed by Plastino and Tsallis 
[7]. This proof-of-principle analysis for generalized VPA (GVPA) models [7] 
established the predominance of concavity of the measure of uncertainty over 
extensivity (as defined within the context of B-G-S statistics). Further work 
by Lenzi, Malacarne, and Mendes [8] and Mendes et. al. [9] demonstrated the 
workings of a GVPA model for the GFE expanded to include second-order 
terms, using a classical harmonic oscillator as an example. More recently, Lu, 
Cai, and Kim [10] demonstrated that the inclusion of higher-order terms can 
significantly improve the results of a GVPA. 

The generic procedure for GVPA's is as follows: (i) evaluation the canonical 
probability distribution pn that maximizes the Tsallis entropy, and, formula- 
tion of the generalized free energy (GFE) [11] 

1 1 - 1 1 

f, = C,-^S, = --^-^ = --ln,Z. (1) 

where, 5*^ is the Tsallis entropy (defined in Section 2), Uq is the generalized 
internal energy, /3 (the energy Lagrange multiplier) is the "inverse thermody- 
namic temperature"!^, and Z is canonical partition function, (u) assuming a 
Hamiltonian of the system as 

H = Ho + XH,, (2) 

where, Hq is a Hamiltonian of a soluble system and XHi is a perturbation to 
Hq. Here, A G [0, 1] is a perturbation parameter, (Hi) perturbation expansion 
of the GFE 

F,(A) = + A5«F, ^ ^ + §6(')Fq ^ ^ + , (3) 



A=0 



A=0 



^ The terms generalized statistics, nonadditive statistics, and nonextensive statis- 
tics are used interchangeably. 

^ A continually updated bibliography of works related to nonextensive statistics 
may be found at http://tsallis.cat.cbpf.br/biblio.htm, 

^ Note that for constraints expressed in the form of normal averages, /3 is replaced 
by /3 = /3/g as is discussed in Section 3.2 and Appendix A of this paper. 
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where: S^'^^Fg — = 1,..., and, (iv) solving (2) with the aid of the 

Hellmann-Feynman theorem [12] 

dfjfl I I 11 /l7'rl\ / a\ 

~dX^\'^\~dX '^^ " ^ '^^ ' 

where, En is the n*'^ energy eigenvalue of the Hamiltonian H, and, {v) evalu- 
ating the generahzed Bogohubov inequality [7-10] 



Fg < F('^ + {H,f\ (5) 

where, {Hi)^^^ is the generalized expectation of Hi in (2) evaluated with A = 0. 
Note that while the proof-of-principle analysis in [7] implicitly assumed the 
above steps (i) — (v) to have been performed, their explicit implementation was 
demonstrated in [8-10]. Further, the expression for the GFE (1) has recently 
been the object of much research and debate. Most generally, the inverse 
temperature /3, which in essence is the Lagrange multiplier associated with 
the internal energy relates to the thermodynamic temperature T as: (3 = -j^, 
where kb is the Boltzmann constant (sometimes set to unity for the sake of 
convenience) only in the limiting case g — ?> 1. Prominent attempts to mitigate 
this issue are those by Abe et. al. [13], Abe [14], amongst others. 

Generalized statistics and the problem of obtaining maximum Tsallis entropy 
canonical probability distributions has been associated with many forms of 
constraints. These include the linear constraints originally employed by Tsallis 
[1] (also known as normal averages) of the form: {A) = Y^PiAi, the Curado- 

i 

Tsallis (C-T) constraints [15] of the form: (A)^ — J^Pi^i i and, the normahzed 

i 

Tsallis- Mendes-Plastino (T-M-P) constraints [11] (also known as g-averages) 

q 

of the form: {{A))^ — J2 ^^^^ ■ The normalized T-M-P constraints render 

i 

the canonical probability distribution to be self-referential, owing to the de- 
pendence of the expectation value on the normalized pdf. A fourth form of 
constraints are the optimal Lagrange multiplier (OLM) constraints [16, 17]. 
A work by Ferri, Martinez, and Plastino [18] introduces a methodology to 
"rescue" the normal averages constraints, and, seamlessly relates canonical 
probability distributions obtained using the normal averages, C-T, ^-averages, 
and OLM constraints. 

The normal averages constraints [1] were initially abandoned because of diffi- 
culties encountered in obtaining an acceptable form for the partition function, 
and, the nonextensive statistics community have largely utilized ^-averages 
constraints instead of the C-T constraints, since (1)^ ^ 1. Recent studies 
by Abe [19,20] suggest that unlike g-averages, normal averages are physi- 
cal and consistent with both the generalized H-theorem and the generalized 
Stosszahlansatz (molecular chaos hypothesis). 
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Despite this physics-related deficiency prominently cited prior works on GVPA's 
[8-10] utilized the C-T constraints as a consequence of mathematical neces- 
sity. For C-T constraints, the generalized internal energy: Ug =< En >q~ 
J2Pn^n ■ The canonical distribution that maximizes the TsaUis entropy is 

n 

1 

[8-10,15]: pn = p{En) = '^"^^"iff"^^ - The partition function is: = 
J] [1 — (1 — g) PEn]'^, where, P is the Lagrange multipher associated with 

n 

the generalized internal energy Uq. 

Prom (1), the GFE is: Eg = -^^^^j^^- Thus, the first-order Newtonian 
derivative of the GFE is of the form: § = ^ (13)"^ E [1 - (1 - ?) PK]^ 

n 

dF 

The first-order Newtonian derivative of the GFE naturally yields: = 
^Pn^lx = ■ The mathematical expression for the q- deformed exponen- 

tial: exp^a; = [1 -|- (1 — g) a;] [21] which, satisfies the differential equation: 
^ = 1/^, mandates that Newtonian derivatives of the GFE (1) results in ex- 
pectations defined in the C-T form. These strictures severely constrain the 
scope and generahty of the GVPA analyses described in [8-10]. 

The primary leitmotif of this paper is to derive a candidate q- deformed GVPA 
employing the q-deformed derivatives defined by [21]: D'f^-^y = [1 + (1 ^ Q") x\ 
where the q-deformed exponential satisfies:D^^^|/ = y. One of the significant 
consequences of such a generalization is the ability to define expectations de- 
fined in terms of normal averages, when taking derivatives of the GFE. 

It has been suggested [10] that the results of Ref. [18] could be employed 
to transform GVPA's derived using C-T constraints to equivalent forms de- 
scribed by either normal averages or g-averages. While such a suggestion may 
be true in principle for simple cases, its practical tractability is questionable 
when applied to studies in communication theory and machine learning, where 
the variational cxtrcmization of the maximum Tsallis entropy problem is done 
with respect to conditional (transition) probabilities. Furthermore, this pre- 
scription does not result in a true generalization of the variational procedure 
employed in GVPA's, within the framework of q-deformed calculus [21]. 

In accordance with [18] , a generic form of the canonical probability distribution 
that maximizes the Tsallis entropy is 



Pn Z{f}*) ' ^ Pn 



(6) 
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For normal averages constraints 



q-' = 2-q 

= ^ = ^ , ; where : /3 = ^L, (7) 

'^q = 'Epl,and,Ug = T,PnEn- 

n n 

For the C-T constraints 

q* = q 

13* = fi, (8) 

n 

For normal averages constraints, (6) and (7) employ the additive duality of gen- 
eralized statistics (see [1] and the references therein). Specifically, g — )■ 2 — g = 
g*IZ]Here, (6) is subjected to the Tsallis cut-off condition [1]: [1 - (f - q*) l3*E„] 
0. The rationale for denoting the generalized internal energy as: Uq, is to high- 
light its nonextensive nature irrespective of the form of expectation involved 
in its definition. Note that even for expectations defined by normal averages, 
Uq has an implicit g-dependency facilitated by the canonical probability (6). 

While replacing the Shannon entropy with the Tsallis entropy yielding GVPA's 
[8-10] represents an initial level of generalization, the q-deformed GVPA model 
described in this paper achieves a further level of generalization by replacing 
the Newtonian derivative with the q-deformed derivative in the variational 
procedure. Specifically, given a function: F (r) = J2F{Tn), the chain rule 

n 

yields: = • Replacing the Newtonian derivative: ^^^7^ by the q- 

deformed derivative defined by [21]: DJ^^F(r) = [1 + (1 — g) r] (see Sec- 
tion 2)and defining: DJ^-^F (r) ^ = S(q)^rF{T), facilitates the transformation: 

^^X^ — )■ 6(q)^rF (t). Note that the q-deformed derivative DJ^-^ operates only on 
the term F{t), which in this paper is a thermodynamic function expressed as 
a q-deformed exponential. Within the scope of this paper, F{t) comprises the 
GFE. Thus the increasing order to derivatives acquire the form 

5(,),.F(r) = /^[,)F(r)f , 



< 



^\:)\Fir) = Dlq^fA,F{r)-k = l,.... 

This results in q-deformed GVPA models that simultaneously exhibit gener- 
alization within the context of the generalized statistics framework employed, 

^ Note that q ^ q* denotes re-parameterization from qio q* . 
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and, in their underlying mathematical structure. The commutation relations: 
Dl,)F (r) = E Dl;^F (rj, and, (r) = DJ^^F (r) f = E Dl^^F (rj ^ = 

J2S(^q-j^T„F (th) are established in Theorem 1. These relations are critical to the 

n 

correctness and admissability of the q-deformed GVPA model, and are derived 
in Section 2 of this paper. 

Section 3 derives a q-deformed GVPA model employing the q-deformed deriva- 
tive in the variational procedure. This is accomplished with the aid of the 
Hellmann-Feynman theorem [12]. Section 3 also analyzes the commutation be- 
tween the Newtonian derivative ^ and the summation sign, and the contribu- 
tion of the cut-off in the variational and perturbation methods for higher-order 
perturbation terms. It is noteworthy to mention that expectations in terms of 
normal averages are achieved by relating the "inverse thermodynamic tempera- 
ture" (3 (the energy Lagrange multiplier) to the "physical inverse temperature" : 
P* [13] via (7), as is demonstrated in Section 3. This is accomplished through 
the relation in (7): /3 = /3*Z{(3*y*~^, where: (3 = (3/q is the scaled "inverse 
thermodynamic temperature" . 

Section 4 formulates the generalized Bogoliubov inequality [5,7-10] truncated 
at first-order terms for the case of the classical harmonic oscillator. Numerical 
examples demonstrating the results for the q-deformed GVPA are presented in 
Section 5. The ability of the q-deformed GVPA model presented in this paper to 
demonstrate both sub-extensivity (sub-additivity) and super- extensivity (super- 
additivity) within the context of the generalized Bogoliubov inequality truncated 
at first- order terms is demonstrated. This feature is not possessed by existing 
GVPA models [8, 10]. Section 6 concludes this paper. 



2 Theoretical preliminaries 



2.1 Tsallis entropy 



The q-deformed logarithm and exponential are defined as [21] 



lug (x) = 
and, 

exp (x) 



l-q ' 



[l + {l-q) x]—^;l + {l-q)x>0 
0; othewise, 



(10) 
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respectively. By definition, the un-normalized Tsallis entropy, is defined in 
terms of discrete variables as [1, 2] 

1 - Ep^ 

S, = = -EPn Pn;EPn = 1- (H) 

? n n 

The constant q is referred to as the nonadditivity parameter. 

2.2 Results from q-algebra and q-calculus 



The q- deformed addition (Bq and the q- deformed subtraction are defined 
as 

x®qy = X + y + {1 - q)xy, ^^^^ 
OqV = i^(=^; l + {l-q)y>0^xeqy^ 

The q-deformed derivative, is defined as 

2.3 Deformed calculus framework for q-deformed GVPA 

Consider a function: F {g (E)) = ^ [1 - (1 - g) (K)]^ = E (£^n))- 

n n 

On comparison with (6), F{t) may be construed as a generic form of the 
canonical partition function Z{f3). Substituting —(3g{E) — r, yields F {t) — 

E[l + (l-g)T„]^=Ei^(r„). 

n n 

1 

Theorem 1: Given a function F (r) = X] [1 + (1 — ^n] ^""^ = E (t™), where 

n n 

T = {ti, , Tat} arc N separate instances of r„; n = 1, A^. The following 

relation involving action of the q-deformed derivative 



Dlq)F (r) = E [1 + (1 - ^n]^ = E i^n), (W) 

n n 

and, 

^ = 5 [1 + (1 - ^ ^ = 5 i^n) (15) 

hold true. 
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Proof: Employing (13) yields 



dFjr) 
dr 



Dl^^F (r) = [1 + (1 - g) r] 

= [l + (l-?)ri]^ + ... + [l + (l-g)r^]^ 
= F (n) + ... + Dl-F M = E (r„). 



(16) 



Similarly, 



D 



+ [l + (l-?)riv]^^ 



(17) 



_d_ 

dX 



dTn 

dX ■ 



Note that (17) assumes that the Newtonian derivative ^ commutes with the 
summation sign [8]. The results of Theorem 1 form the basis for the q- deformed 
GVPA model presented in this paper, which extends earlier GVPA studies [8- 
10]. Note that (17) may be extended to obtain expressions for higher order 
derivatives with the aid of (9), which is demonstrated in Section 3. Setting: 
F (r) = [1 + (1 - g) r]^ = E [1 + (1 - q) Tn]^, DT.F^t) yields with the 
aid of (16) 

Dl,)F (r) = [1 + (1 - g) r] ^ 

= E [1 + (1 - g) Tn] [1 + (1 - q) Tn]^ (18) 

n 

= E[l + (l-g)rn]^ =Ei^(rn). 



Analogously 



dF (r) 
dr 



= 5:[l + (l-g)T„]^. 



(19) 



Setting r^-^E^%^ <5(,),.F (r) = Dl^^F (r) g = Dl^^F (r) §f = 
-pDl^^Fir) ^. Thus, (17) and (18) yield 



= -/3E[i-(i-g)/3K]^^ 



A=0 



(20) 



A=0 



Analogously 



dF (t) 



dX 



x=o 



(21) 



A=0 



8 



Comparison between (20) and (21) provides initial evidence that optimality 
conditions obtained using the q- deformed derivative quahtatively differ from 
those obtained using the Newtonian derivative. 



3 q-Deformed generalized variational perturbation approximation 



3.1 Overview of the Rayleigh-Schrddinger perturbation method 



Consider the eigenvalue equation 



H\n) 



En\n) 



(22) 



where H is the Hamiltonian, En is the energy eigenvalue, n is the quan- 
tum number, and, \n) is the eigenfunction (eigenvector). Note that the non- 
degenerate case is considered herein. The eigenf unctions and the eigenvalues 
are written as [22] 



\n)^\nf^+X\nf^ + X-'\nf^ + ..., 
En = E^) + \EW + X^E^) + ..., 



(23) 



where A G [0, 1]. In (23), E^^^ and \nY^' are the k*^ corrections to the energy 
eigenvalues and eigenfunctions, respectively. Note that the non-degenerate the- 



Substituting (23) 



ory IS valid only when: K") (m| Hi |n)^°^| « - 
into the time-independent Schrodinger equation, and comparing the coeffi- 
cients of the powers of A on either side, the first two eigenfunction corrections 
are 



\n 



(2) 



^ (0)(m|gi|n)W I UP) 



,(0) 



E 



Hi Hi 

nn 



1 (-^Imn) 



\m) 



(0) 



(24) 
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3.2 Generalized free energy correction terms 



Setting T„ = —P*En, (6) acquires the form 



Pn = 



[i+(i-g*)r„]i- 
Z{r) 



(25) 



For expectations defined by normal averages, tfie GFE (1) is defined by: 
Fq — Uq — ^Sq — —^luq Z (^^^ , whcrc: P — P/q. Note tfiat the scaled energy 

Lagrange multiplier /3 is introduced in order to achieve consistency between 
the two forms of the GFE. The q* -deformed GFE is defined as 



Oq* 



(26) 



P q*-l 

The derivation of (26) is provided in Appendix A of this paper. Defining 

dr 

Si,^),rZ {t) = Dlq,^Z {t) -, (27) 
a perturbation expansion of the q* — deformed GFE, analogous to (3), is 



Oq* — 'J{q*),T^q* } 



A=0 



C-{2) 
2! ^q* 



A=0 



+ ... 



(1) 



(28) 



3. 3 Derivation of terms in the GFE perturbation expansion 



From (23) and (26), the zeroth-order term is 



(0) ^ _ 1 ^(,*),o(r)"*-^-l 
9* ^ q*-l ■ 



where 

V),o(^*) = E[i-(l-g*)^*£^f 



1 



(29) 



Given the canonical probability distribution and the canonical partition func- 
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tion in (25), (28) yields 



= [l + (l-,*)r]^^g 



(30) 



Substituting (25) and (26) into (30), and employing (15), yields 

^-IZ irf-' x: [1 + (1 - g*) r„] 4^ [1 + (1 - g*) r„] ^ ^ 
-i^(r)^*-'E[l + (l-g*)r„]T^^ 



Setting: r„ = -j3*En, employing the chain rule: ^ = substitut- 
ing: ^ = P*Z{P*Y-^ from (7) into (31) yields 

(32) 



Employing the Hellmann-Feynman theorem (4), the first-order perturbation 
term in (28) is defined as 



(1) 



A=0 



^E[i-(i-g*)/S*K]^^ 



A=0 



^Efi-(i-?*)r£^f 



A=0 



77. ^ ' 



(0) 



(33) 



where, {Hi)^^J denotes the expectation of Hi in normal averages form, defined 
with respect to the probability p{E^^) parameterized by q*. 

Note that $>qVin (32) is obtained by operating on the q* -deformed GFE in (26) 
with S(q*)^T-n, where r„ = —j3*En. Thus, the "inverse thermodynamic temper- 
ature" /3 (the energy Lagrange multiplier) is related to the "physical inverse 
temperature" /3* by the relation: = /3 = P*Z{(5*Y employing the re- 
sults of Theorem 1. In (33), $5^^ ^ is tantamount to specifying the En = E^^ 

as defined in (4) and (23). Thus, /3 7^ Z^q*)^Q{P*Y'-^ , where Z(q.),o(/9*) is de- 
fined in (29). 
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Setting Tn — —(3* En, (32) is re-expressed as 
Operating on (34) with 5(q*),r defined in (28), yields 

C>(2) _ ^ c>(l) _ T^T c>(l) dr 

Employing the Leibnitz rule for deformed derivatives [21] 

Dl^.^ [A (r) B (r)] = B (r) Dl^.^A (r) + A (r) DJ^.^B (r) , 

and defining: A (r) = E [1 + (1 - ?*) ^n]^ ^ and B (r) = Z (t)"\ 

^ V ' 

Term~l 

~{i:ii+(i-«>»i'^^^[,-,^M-')|. 

^ V ' 

rer-m-2 

Evaluating (37) term-wise with the aid of (15), yields 

Term - 1 = -X |^ (^)-i ^r^^^ E [1 + (1 - g*) Tn]^^} g 
= (r)-^ E [1 + (1 - q*) Tn]^ (^)' 

+ Z (r)"^ E [1 + (1 - r) r.] ^ 4: (^) (^) } 

= (^)"' 5 [1 + (1 - ^r^]^ 

+ Z (r)-' E [1 + (1 - q*) Tnf^^A (l^) (1^) } 
= -j^ [Z {rV E [1 + (1 - r) Tn]^ (^)' 
+ ^(r)-^E[l + (l-r)r.]^^j, 



and 
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Term - 2 = |e [1 + (1 - g*) r„] ^ ^ (r)"^] | g 

= ^ (r)-^ E [1 + (1 - g*) r.]^ (^)'} . 



2^ ^2, 



+^(r)-^E[l + (l-g*)r„]T-?^ 

+ ^ (r)-' E [1 + (1 - r) Tn]^ (^)'} 



(39) 



(40) 



Setting Tn = -/3*^n, (40) yields 



c>(2) 



-Z {(3*)-' 5 [1 - (1 - (ff )'} 



+z (/3*)-^ E [1 - (1 - g*) 



(41) 



With the aid of the Hellmann-Peynman theorem (4), the second-order pertur- 
bation term in (28) is 



C>(2) 



A=0 



r). ^ ' 



(42) 



A=0 



Employing (23) ^_^= \n)^^^) and (24), (42) acquires the compact form 



C>(2) 



A=0 



(43) 



where: i?i,„„ ='•''■' (n| ifi |m)*'*^''. Note that (43) is the direct equivalent of Eq. 
(6) in Ref. [8] and Eq. (24) in Ref. [10]. In the classical limit, the last term 
in (43) vanishes, since the Hamiltonians and Hi commute. In keeping 
with the gist of previously cited works on GVPA's [8,10], this paper employs 
the classical harmonic oscillator to examine the properties and efficacy of 
the q-deformed GVPA principle introduced herein. Thus, the second-order 
perturbation term becomes 



C>(2) 



A=0 



(44) 
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Thus, for commuting Hamiltonians Ho and Hi, the perturbation expansion of 
the GFE is 



v(A)=- f"-"rr''' +A(j/.);°' 

+§ir{{{H,)f]y--ZPn{Eif){H,,j] + 0{X') + ... 

The above analysis and (17) make tacit assumptions concerning the com- 
mutation between the Newtonian derivative ^ and the summation sign. To 
analyze the contribution of the cut-off in the q-deformed GVPA model for 
P* > and En > 0, parallels are drawn with the analysis in Ref. [8]. First, 
the summation sign J2 is replaced by the integral / dV, where: F = H 

n s " 

Ps is the canonical momentum, and h is Planck's constant. The integration 
is performed over the phase-space region defined by: [1 — (1 — q*) /3*H] > . 
Defining: / = [1 - (1 - q*) f3*H]^ = [1 + (1 - q*) = /(r), the follow- 

ing identity is employed by invoking (9) and Theorem 1 

i,IdrDl^.)f{r)^JdrDl^.^f{r)§, 
+ Ievi:dS^{Dl;,^fir^)^/i-). 



Note that the distinction of (46) vis-d-vis Eq. (10) in Ref. [8] is due to re- 
placing: by: 5(q*),T/(T) (see (9)), in the q-deformed GVPA model. Here, 
/(r) is a function of the phase-space variables and, dV is the hypersur- 
face defined by: [1 - (1 - q*) P*H] = [1 + (1 - g*) r] = 0. The n^'^-order per- 

/ \ n— (n — l)g 

turbation 0=^. contains a term proportional to: [H- (1 — g*)r] ^ . For 

o * 

2 — q 

example, Qg*' in (40), contains a term: [1 -|- (1 — g*) t„] corresponding 

n 

to n = 2. In order that 9^*'* not contribute to the second term in (46), the 
condition: n — {n — 1) q* > q* < 1 H — ^ is to be observed so that: 
Dl,.)f (r) = ^ / (r) = on dV. 



4 The generalized Bogoliubov inequality 



Within the context of the q-deformed GVPA, the generalized Bogoliubov in- 
equality truncated to first-order terms is 

< + {Hi)P . (47) 



For a 1-D classical harmonic oscillator of mass M and angular frequency u, 
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the Hamiltonian is: H = ^ + Mjiij^_ Here, x is the coordinate and p is the 
canonical momentum. Following a procedure analogous to that in Ref. [10], 
from (6), the canonical partition function is expressed in continuous form as 

JV 



Z in = / [1 - (1 - P*n5o]^dn = P* > 0. (4J 







Here, < g* < 2, oo for g* > 1, and = for q* < 1. Akin 

to [10], 6o = fkj, is a positive constant with units of energy. Also, h = /i/27r, 
where h is Planck's constant. Thus, (26) is re-written as 



3„ 



(49) 



The unperturbed Hamiltonian for a particle of mass M in a 1-D box is: i^o = 

2 

2M+Vo, where Vq = for |x| < L/2 and Vq — )■ oo for |x| > L/2. In accordance 



with [10], the continuous energy spectrum of the particle is: — 2ml2, 

where 6^ = hn is a constant with the dimension of action. Thus, (29) is 
described in terms of the Euler F function as 



N 







Mt, 



{l-q*)P*^ 



2ML2 



2{l-q*)p* r(f-r — ^ 
L I Mtt ' '"(~i + ?^ 

2(r-i)/3* r(^) 



1 \ ) y ^ 



1 



dn 



-;3 > g* > 1 



(50) 



where A^ ^ cx) for g* > 1 and A^ = y\Ijtz^ for g* < 10 Note that (50) is 



identical to Eq. (11) in [8], with g* replacing g. 



Following the procedure employed in [10], the matrix elements of in (33) 
and (43) are (Eq. (33) in [10]) 



L/2 



Hi. 



^n.m2 



-L/2 



Mco X dx 
2 2L 



■^2A 



(51) 



Here, 5nm is the Kronecker delta. Employing (50), the first-order perturbation 



Symbolic integration was performed using MATHEMATICA® . 
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term (33) becomes 



c>(l) 



A=0 



24 



N 



Hi*) 



{l-q*)P' 



v?5l 
2ML2 



1 



dn 



2t2 



24 



(52) 



From (25), (26), (29), (48)-(50), and (52), the generalized Bogoliubov inequal- 
ity (47) is 



1 



Z(/3 



- 1 



1 



(g*-l)/3 



U2-9*)/3* 



2. \ 



< ^ 

where 



Z{q*),0 (/?*) 



+ 



2 7-2 



24 ' 



L / 2Af 



V),or)= ^V(i-.-)rr(F;X^' 

L I 217;? U~2 + ?^j 



ii'^^W r(^) 



3 > g* > 1 



and 



(53) 



The value of L is obtained by minimizing the right hand side of (53), yielding 



/3 ^1 

where 



12 



^ L 



12C„.(/3*) 



1 

9^ 



2Mtt 



;g*<l 



(54) 



> . 



2Mtt 



2 ' q* -1 , 



; 3 > g* > 1 



Here, (54) demonstrates that L accounts for both suh-additivity [q* > 1) and 
super- additivity [q* < 1). This feature is absent in formulations of the gen- 
eralized Bogoliubov inequality truncated at first-order terms in Refs. [8, 10] 
Specifically, Eq. (13) in Ref. [8] and Eq. (40) in Ref. [10] make no allowance 
for sub- additivity and super- additivity in the expression for L. 



Note that sub- additivity and super- additivity are defined in terms of q* because 
(53) and (54) are parameterized by q*. 
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5 Numerical studies 



The q-deformed GVPA model is numerically studied by evaluating the gen- 
eralized Bogoliubov inequality (53) for values of /3* e [0.01,3.5], akin to 

the parametric perspective described in [18] for various values of q*. Here, 
M = u = h = 1. Representative examples for the generalized Bogoliubov in- 
equality (47) and (53) are demonstrated in Figure 1, for q* = 0.5 and q* = 0.95. 
Numerical examples for the generalized Bogoliubov inequality for g'* = 1.3 and 
q* — 1.75 are depicted in Figure 2. 

From Figure 1 and Figure 2, it is readily observed that the generalized Bogoli- 
ubov inequality in (47) and (53) tends to an equality with decreasing values of 
q*. Specifically, the difference between the exact solution {^q*) and the pertur- 
bation solution (3^^?+ < Hi >f*) decreases with decreasing q*. Note that the 
exact solution of the GFE (the LHS of the inequality in (47) and (53)) almost 
exactly coincides with the perturbation solution, for q* = 0.5. The relations 
between (3 and (3* are displayed in Figure 3 for both q* < 1 and q* > 1 . 

The expression for L (54) is the most explicit manifestation of sub-additivity 
and super- additivity possessed by the q-deformed GVPA, for k^^-oidei pertur- 
bation expansions of the q* — deformed GFE (26), k > 1. Figure 4 displays 
the dependence of L on For < g* < 1, L increases with increasing values 
of q* . In contrast, for 1 < g* < 3, L decreases with increasing values of q* . 



6 Summary and conclusions 

The theoretical framework for a q-deformed GVPA model which generalizes 
previous works on GVPA models [8-10], has been formulated. The underlying 

theory for the variational procedure of the q-deformed GVPA model employs 
q-deformed calculus [21]. The significant feature of the q-deformed GVPA is 
that expectation values may be self-consistently formulated in the form of nor- 
mal averages, instead of the Curado-TsaUis form [8-10]. This feature acquires 
special significance owing to recent results that establish that expectations in 
the normal averages form, in contrast to the vastly more utilized g-averages 
form, are physical and consistent with both the generalized H-theorem and 
the generalized Stosszahlansatz (molecular chaos hypothesis) [19, 20]. 

It is qualitatively and quantitatively demonstrated that the q-deformed GVPA 
model exhibits both sub-additivity and super- additivity in terms of the non- 
additivity parameter q* for the generalized Bogoliubov inequality, truncated 
at first-order terms. This property is not possessed by previous GVPA models 
[8, 10]. Specifically, it may be construed that the q-deformed GVPA presented 
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in this paper exhibits authentic characteristics of a VPA even for the general- 
ized Bogohubov inequahty truncated at first-order terms. In contrast, previous 
cited analyses [8, 10] do not exhibit any equivalent property, rendering them 
purely variational principles when the generalized Bogoliubov inequality is 
truncated at first-order terms. Numerical simulations that demonstrate the 
results and the efficacy of the q-deformed GVPA are presented. 

Future works that will be presented elsewhere accomplish a three-fold ob- 
jective: (i) a comparative study taking into account higher-order terms be- 
tween the q-deformed GVPA model presented in this paper, its counterpart 
based on the dual Tsallis entropy: Sq*=2-q (see Ref. [3] and the references 
therein), and the results of previous cited studies [8-10], (ii) a q-deformed 
GVPA model for the homogeneous Arimoto entropy [23] with normal averages 
constraints, defined in terms of the escort probability [24]. The homogeneous 
nonadditive Arimoto entropy is defined in terms of the escort probability as: 

(P) = — ^ ' , and, (Hi) use of the q-deformed GVPA to analyze 

critical point behavior in deterministic annealing [25], within the framework 
of generalized statistics, and, the deformed statistics information bottleneck 
method [4] in machine learning. 
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FIGURE CAPTIONS 



Fig. 1: Generalized Bogoliubov inequality for the q — deformed GVPA model 
for q* = 0.5 and q* — 0.95. Here, Qg* is the LHS of (53) (the exact solution) 
and $5^?+ < Hi >fJ is the RHS of (53) (the perturbation solution). The 
generalized Bogoliubov inequality increases with increasing q*. Note the near 
overlap of the exact solution and the perturbation solution for q* = 0.5. 
Fig. 2: Generahzed Bogoliubov inequahty for the q — deformed GVPA model 
for = 1.3 and q* = 1.75. Here, ^g* is the LHS of (53) (the exact solution) 
and + < Hi is the RHS of (53) (the perturbation solution). Note the 
increase in the generalized Bogoliubov inequality with increasing q*. 
Fig. 3: Relation between the scaled "inverse thermodynamic temperature" 
{P/q = P) ^-iid the "physical inverse temperature" (/?*) for g*=0.5, 0.95, 1.3, 
and, 1.75. 

Fig. 4: Dependence of L on Note the increase in the value of L for in- 
creasing q* <1, and the decrease in the value of L for increasing q* > 1. 
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Appendix A: Derivation of expression for gf*-deformed Generalized 
Free Energy 



The canonical probability that maximizes the Tsallis entropy using constraints 
defined in terms of normal averages is (Eq. (23) in Ref. [18]) 



Pn = 

where 



(A.l) 



Ug = EPnEn, P=l,^g = EpI and, Z (^) = (k, + {q - 1) pUg) 



Prom the definition of the partition function in (A.l), the following thermo- 
dynamic relation is obtained 



' l-q dp 



= -U,. (A.2) 



The Tsallis entropy is 



l-q 



(A.3) 



Using the definition of Z{P) in (A.l), (A.3) yields the thermodynamic relation 

1-9 



+ (l-q)(3U,-l dS, 



l-q 



dU, 



' -p. 



(A.4) 



Prom (A.4), the GPE is thus defined as 



Fq — Uq — ^Sq 



Z{B) '+(l-g)/3f/,-l 
l-q 



-iln,z(/3). 



(A.5) 



Prom (A.l) and (A.5) 



l-{q-l)j^.E^ 



1 

9-1 



-i In, E [1 - (g - 1) P*E^]—^ = -i In, Z iP*) , 



(A.6) 
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where: (3* - Here, Z {/3*) = E [1 - (? - 1) 

Invoking the additive duality: q* = 2 — q (see Ref. [1] and the references 
therein) resulting in: Z{p*) = E [1 - (1 - q*) P*Er^]^, setting: t„ = -p*E^ 

n 

and employing Eq. (25) of this paper, the q* -deformed GFE (Eq. (26) of this 
paper) is expressed as 

i^.-.,*=2-, - - -J - -J _ ^ . (A.7) 



26 



